Feedback and Rate Asymmetry of the Josephson Junction Noise Detector 



OO 
O 

o 

(N 
•*-> 

o 
O 



c 

o 
o 



> 

OO 

m 

On 
(N 

d 

00 

o 



X 



D. F. Urban and Hermann Grabert 
Physikalisches Institut and Freiburg Institute for Advanced Studies, 
Albert-Ludwigs-Universitat, 79104 Freiburg, Germany 
(Dated: October 16, 2008) 

The Josephson junction noise detector measures the skewness of non-Gaussian noise via the asym- 
metry of the rate of escape from the zero- voltage state upon reversal of the bias current. The feed- 
back of this detector on the noise generating device is investigated in detail. Concise predictions 
are made for a second Josephson junction as noise generating device. The strong nonlinearity of 
this component implies particularly strong feedback effects, including a change of sign of the rate 
asymmetry as the applied voltage approaches twice the superconducting gap. 

PACS numbers: 85.25.Cp, 72.70.+m, 73.23.-b 



In the last decade there have been extensive theo- 
retical studies of non-Gaussian noise generated by non- 
linear electronic nanostructures^. This new field of 
full counting statistics (FCS) has put forward numer- 
ous predictions indicating that the higher order noise 
cumulants contain valuable information about the elec- 
tronic transport mechanisms within the noise generat- 
ing element that are not accessible from more standard 
measurements of the conductance and the noise power 
alone. A few years ago, first experimental observations 
of non-Gaussian current noise have been reported 3 -^, and 
more recently quantum point contacts 5,6 and Josephson 
junctions^ have emerged as promising on-chip noise de- 
tectors. 

The current experimental efforts mainly focus on re- 
liable measurements of the third noise cumulant C3, the 
skewness of the noise. A Josephson junction (JJ) in the 
zero- voltage state can detect the skewness of a noise cur- 
rent passing the detector since C3 is odd under time rever- 
sal and thus leads to an asymmetry of the rate of escape 
from the zero- voltage state of the JJ when the direction of 
the bias current is reversed^ After some primary sug- 
gestions for JJ noise detector o 9 ' 10 ! 11 ! 12 ! 13 , concrete the- 
oretical studies for the setup used in experiments have 
started recently 14 i 15 ' 16 

Now, that the theoretical tools for predictions for JJ 
noise detectors with realistic parameters are available, 
we shall address here the feedback of the noise detector 
on the noise generating device and show that this feed- 
back can modify the rate asymmetry, the quantity deter- 
mined experimentally quite considerably. A detailed un- 
derstanding of this feedback is essential for reliable data 
analysis. But this issue is also of relevance for all on-chip 
detectors, since this very concept implies two mesoscopic 
devices, detector and measured device, that interact on 
the same chip and thus have to be treated on the same 
footing. To be concrete, we shall specifically investigate 
the detection of non-Gaussian noise generated by a an- 
other JJ biased by a voltage Vat with eV/v close to twice 
the gap A of the superconductor. Because of the strong 
nonlinearity of the noise generating device in this region, 
feedback of the detector leads to particularly pronounced 
modifications of the rate asymmetry. 
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FIG. 1: Circuit diagram of a JJ noise detector: A JJ with 
critical current I c and capacitance C is biased in a twofold 
way. The branch to the right puts an Ohmic resistor Rb 
is series with the junction and is biased by the voltage Vs. 
The branch to the left is biased by a voltage Vn and includes 
the noise generating nonlinear element characterized by its 
cumulants C,-. 



A typical realization of a JJ noise detector is sketched 
in Fig. [1] A JJ with capacitance C and critical current I c 
is connected by two branches which act as noise sources. 
A bias voltage Vb is applied to the first branch with 
an Ohmic resistor Rb in series with the junction. This 
branch allows to control the detector by fixing the av- 
erage bias current, applying current pulses, and reading 
out voltage signals. A second voltage Vjv is applied to 
another branch which includes a non-linear noise gener- 
ating conductor, again in series with the detector. This 
noise generating element can most generally be charac- 
terized by its charge transfer statistics, respectively its 
voltage dependent current cumulants {Ci(V),C2(V), ...}. 
Current experimental set-ups are more sophisticated, but 
the circuit diagram in Fig. [1] captures the essentials of a 
JJ on-chip noise detector. 

The state variables of the J J detector are the charge Q 
on the junction capacitance C and the phase difference cp 
between the order parameters of the superconductors on 
either side of the tunnel barrier. The phase dynamics of 
the detector can be described by classical physics, since 
the effective temperature determined by the strength of 
Gaussian noise from the two branches (see below) is typ- 
ically much larger that the crossover temperature below 
which macroscopic quantum tunneling of the phased be- 
comes relevant. The stochastic dynamics of the setup can 
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then be described in terms of a Hamiltonian H— which 
depends on the state variables Q, ip and the conjugate 
thermodynamic forces A, [i. We note that the Hamil- 
tonian of the stochastic theory has to be distinguished 
from the Hamiltonian of the microscopic theory. TL de- 
pends only on gross variables of the system and is mea- 
sured in the same units as the rate of change of entropy 
[fcgu] rather than energy [huj]. For the setup in Fig. [1] 
H = Hd + 7~Cb, where 

2e Q 

'Hd = y^^i - / c sin(</j)A (1) 
describes plasma oscillations of the JJ detector, and 



includes the leading order effects of non-Gaussian noise. 
Here we have defined the total bias current through the 
detector 
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and the parallel (differential) resistance 
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of the circuit. Moreover, we have introduced the effective 
temperature 
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adds the two biasing branches. While the branch 
with Ohmic resistor Rb generates Johnson-Nyquist 
noise and has only two nonvanishing cumulants C-f = 
{V B - Q/C) /R B and Cf = 2k B T/R B , the other branch 
is characterized by the cumulants Cj that are taken at 
voltage Vn — Q/C. 

The canonical equations of motion following from the 
Hamiltonian Tt determine the most probable path be- 
tween a given initial state i and a final state /, and the 
transition probability between these states may be writ- 
ten as a path integral^ 



Pi^f = J D[Q,(p, A,/i]exp|-^-A[Q,^,A,/i]| 
where the action functional is given by 



(3) 



A[Q,cp,X,fx] 



ds 



Q\ + <pfj,-H(Q,ip,\,n) . (4) 



During escape of the detector from the zero-voltage 
state, the dimensionless quantity eX/ks can be shown-^ 
to be always much smaller than 1. Furthermore, the 
voltage across the JJ detector 
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that builds up under noise activation is proportional to A. 
Keeping in Eq. @ only terms that are at most of third 
order in A and Q/C, we obtain H = 7Y 2 + H3, where 

H 2 = — -fi-I c sm((p)X+[ I t - — — ) A+^-A (6) 
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describes the stochastic dynamics in presence of Gaussian 
noise, while 
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characterizing the strength of Gaussian noise. The 
are the cumulants Cj taken at voltage Vat. 

In terms of the effective parameters introduced above, 
the leading order Hamiltonian TL2 has the standard form 
for a biased JJ described by the resistively and capaci- 
tively shunted junction model^ The phase cp moves in 
the "tilted washboard" potential 
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[cos(^) + Sip] 
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driven by Gaussian noise. Here s = It/I c is the dimen- 
sionless bias current. For < s < 1, the potential has 
extrema in the phase interval [0, 27r] at 



¥>wcii,to P = arcsin(s) = - =F $ , 
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with S « y/2(l - s) for 1 - s < 1. When the JJ is 
trapped in the state y wc ii = f — S, the average voltage 
Vj across the junction vanishes. This zero-voltage state 
is metastable, and to escape from the well, the junction 
needs to be activated to the barrier top at ip top = § + S 
by noise forces. The weak non-Gaussian noise described 
by H3 also gives a contribution to this process, and the 
rate of escape T may be written as 



r = /e 



-(B2+B3) 



(13) 



where / is the prefactor of the rate, and the expo- 
nential factor is determined by the action of the most 
probable escape path, which is a solution of the canoni- 
cal equations of motion following from the Hamiltonian 
with the boundary conditions (p(t = —00) = ip wc i\ and 
ip(t = +00) = (ftop- The exponential rate factor has a 
dominant contribution 
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arising from TL2 , where 



AU = U(ipt op ) - U(ip we n) 



(14) 



(15) 



3 



is the barrier height of the metastable well. The cor- 
rection B 3 due to non-Gaussian noise may be evaluated 
by treating H3 as a perturbation. Following the lines of 
reasoning of Ref. [H one finds 
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is the effective third noise cumulant. The second and 
third terms in Eq. (|17p arise from the feedback of the JJ 
on the noise generating device, which is a consequence of 
the finite voltage Vj arising during escape. The quantity 
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is expressed in terms of the deterministic trajectory 
Vreiax(i) starting at t = —00 at the barrier top ipt op and 
relaxing (in the absence of noise forces) to the well bot- 
tom tp we \i reached at time t = 00. 

The coefficient B2 is even under time reversal, while 
B3 is odd, like the third noise cumulant. Accordingly, 
when the direction of the bias current is reversed, the 
rate coefficient shows an asymmetry which is a signature 
of non-Gaussian noise. The deviation of the rate ratio 
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= ex P [-2B 3 (I t )]«l-2B 3 (Jt): 



(19) 



from 1 allows for an experimental determination of B3 
despite the fact that the corrections due to non-Gaussian 
noise are typically small. 

In the following we will focus on a specific noise gener- 
ating element, namely a J J with normal state tunneling 
resistance Rt and capacitance CV biased by a voltage 
Vat which is tuned close to twice the gap A of the super- 
conductor. In view of the lead to the voltage source, 
this noise generating element is in series with an ap- 
proximately Ohmic lead resistance Rn in the range of 
a few 100 n. A corresponding setup is currently stud- 
ied experimentally, 2 - and is well suited to examine the 
feedback of the detector on the noise generating ele- 
ment, since the behavior of the JJ is highly nonlinear 
for eV N w 2A. 

A JJ with parameters around those indicated in Tab. 
Uhas a very low sub-gap current, and for eV > 2 A the 
current is carried by tunneling quasiparticles. We can 
then employ the semiconductor modelj^i where the su- 
perconductor is described in terms of quasiparticles with 
a density of states (DOS) N$(E) given by the familiar 
BCS expression 
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A = 200 /ieV superconducting gap 
Cjv = 1 fF capacitance of noise generating JJ 
Rt = 30 kfi tunneling resistance of noise generating JJ 
Rn = 200 Q lead resistance of non-Gaussian noise branch 
Rb = 200 Q lead resistance of biasing control branch 

TABLE I: Circuit parameters used for the calculations of the 
quantities displayed in Figs. [2] and [3] 
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FIG. 2: J-V-characteristic of the noise generating JJ at 
T=50 mK (solid line) and T = 0K (dashed line) in presence 
of DCB. The dotted line shows the step-like discontinuity for 
a JJ without DCB. The inset displays the P(E) function. 



where Ao is the normal state DOS at the Fermi level. 
The mean current through the JJ can be written as 



(/) = 0^ = e(T^ - IV) : 



(21) 



where the forward and backward quasiparticle tunneling 
rates across the junction interface are given 
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and r_(V) = r_>(— V), respectively. Here, f(E) is the 
Fermi distribution. The influence of the electromagnetic 
environment is described in terms of the probability P{E) 
that a tunneling quasiparticle looses the energy E to the 
environment. This dynamical Coulomb blockade (DCB) 
effect leads to a smearing of the expected discontinuity 
of I(V) at eV = 2A. Since the phase dynamics of the 
JJ detector is very slow on the time scale of quasiparticle 
tunneling in the noise generating J J, the impedance of 
the electromagnetic environment is essentially given by 
the resistance Rn of the lead to the voltage source Vn- 
For this case of an Ohmic environment one has22 
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FIG. 3: The effective third cumulant for T=50 mK (solid line) 
and T — K (dashed line). The inset shows the effective 
temperature. 



where Z e g(o;) = i?jv/(l + iojRnCn) is the effective 
impedance and Rk = h/e 2 . For the set of parameters 
given in Tab. HI the resulting /-^-characteristics of the 
noise generating JJ is shown in Fig. [21 The inset shows 
the function P(E). 

The cumulants of the quasiparticlc tunneling current 
are readily determined from the tunneling rates (|22[) . 
One has2^ 

Cf = e J [r_ + (-l) J r^] (24) 

which implies, in particular, Cg P = e 2 C^ p . Although Rn 
has a significant effect on the noise cumulants via DCB, 
in the region 6Vm > 2A, its influence on the voltage 
across the noise generating JJ can safely be neglected 



since R^/Rt <C 1. Hence, the modifications of this volt- 
age are essentially due to the detector feedback of order 
Vj/Vn discussed above. Accordingly, in Eqs. (JTUJ) and 
dlZD , we can put Cf ~ Cf (V N ). 

The highly nonlinear behavior of the cumulants for 
eV ~ 2A strongly affects the effective temperature T e g 
and the effective third cumulant C3.cs. Both quantities 
are shown in Fig. [3] as functions of the applied voltage, 
again for the representative set of parameters given in 
Tab. U While the skewness of the noise C3 is always pos- 
itive for positive It, the rate asymmetry (|19l) . which is 
proportional to C3.cs, passes through zero as Vn is de- 
creased and it takes large negative values for voltages 
slightly above 2A/e. This change of sign constitutes a 
pronounced feedback effect that should be easily observ- 
able experimentally. 

In summary, we have shown that feedback of the JJ 
noise detector on the noise generating device can be very 
pronounced if the device under investigation is highly 
nonlinear. In particular, another JJ as noise source 
should allow for a detailed experimental study of this 
effect, which leads to a strongly enhanced asymmetry of 
the switching rate of the detector JJ and even a change 
of sign of the apparent noise skewness C3 je ff as the volt- 
age applied to the noise generating junction approaches 
twice the superconducting gap. The size of the feedback 
corrections sensitively depends on DCB of quasiparticle 
tunnelling. 
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